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We study static magnetic monopoles in the context of varying-a theories and show that there is 
a group of models for which the t'Hooft-Polyakov solution is still valid. Nevertheless, in general 
static magnetic monopole solutions in varying-a theories depart from the classical t'Hooft-Polyakov 
solution with the electromagnetic energy concentrated inside the core seeding spatial variations of 
the fine structure constant. We show that Equivalence Principle constraints impose tight limits on 
the allowed variations of a induced by magnetic monopoles which confirms the difficulty to generate 
significant large-scale spatial variation of the fine structure constant found in previous works. This 
is true even in the most favorable case where magnetic monopoles are the source for these variations. 



I. INTRODUCTION 

The possibility that at least some of the fundamental 
"constants" of Nature might be dynamical appears nat- 
urally in the context of models with extra-spatial dimen- 
sions. The interest in this type of models has recently 
been increased with results coming from both quasar ab- 
sorption systems 0,0 (see however 0,0,111) and the Oklo 
natural nuclear reactor Q suggesting a cosmological vari- 
ation of the fine-structure constant, a = e 2 / (4 it Tic) at 
low red-shifts. Other constraints at low redshift include 
atomic clocks 6] and meteorites |2J while at high redshifts 
there are also upper limits to the allowed variations of a 
coming from either the Cosmic Microwave Background 
or Big Bang Nucleossynthesis [2l llfj, 111! Il2l \1A 114 . . 

On the theoretical side some effort has been made 
on the construction of self-consistent phenomenological 
models for space-time a variabilit y m ost of them based 
on the orig inal Bekenstein model llal . In some of these 
models [ij, E E HJ El El 0| the possibility that 
the variation of fine structure constant might be coupled 
to the dark energy equation of state responsible for the 
recent acceleration of the Universe [2^, [2(| has also been 
explored. 

Most of these studies have focused mainly on the vari- 
ation of a with time. The reason for this is that it is very 
difficult to generate significant large-scale spatial varia- 
tions of the fine structure constant. This has been shown 
in [27j (see also j^) in a model where the spatial vari- 
ation of a were induced by cosmic strings cosmic and in 
|29j in a more general context. In the cosmic string case, 
the electromagnetic energy concentrated along the core 
of the local string is a source of these spatial variations, 
which are roughly proportional to the gravitational po- 
tential induced by the strings. They are constrained by 
Equivalence Principle to be small and overall limits for 
the allowed spatial variations have been calculated. 

In this article, we consider the case of spatial variations 



of a seeded by a static non-Abelian local monopole. We 
show that the electromagnetic energy in its core can pro- 
duce spatial variations of the fine-structure constant in 
the vicinity of the monopole in the context of Bekenstein- 
type models. We start by solving the standard problem 
with no a variation and recover the classical solution pre- 
sented by t'Hooft and Polyakov HilEl ( see f° r instance 
Refs. IH, 0, HU) and then proceed to investigate 
other solutions with variable a. 

The article is organized as follows. In Sec. II we in- 
troduce Bekenstein-type models in Yang-Mills theories 
and obtain the equations that describe a static magnetic 
monopole. We write the energy density of the monopole 
and show that the standard electromagnetism is recov- 
ered outside the core if a is a constant. In. Sec. Ill 
we present the numerical technique applied to solve the 
equations of motion and give results for several choices of 
the gauge kinetic function. In Sec. IV we study the lim- 
its imposed by the Equivalence Principle on the allowed 
variations of a in the vicinity of the monopole, as a func- 
tion of the symmetry breaking scale. Finally, in Sec. IV, 
we summarize and discuss our results. Throughout this 
paper we shall use units in which h = c = 1 and the 
metric signature H . 



II. BEKENSTEIN-TYPE MODELS IN 
NON-ABELIAN FIELD THEORIES 



We first introduce Bekenstein-type models in a non- 
Abelian Yang-Mills theory. We take the electric charge to 
be a function of the space-time coordinates, e = eo e(cc M ) 
in which e is a real scalar field and eo is an arbitrary 
constant charge. Let the Higgs field <J> a be an isovec- 
tor, where a = 1,2,3 are internal indices associated to 
isospace with SU{2) symmetry. 

The Lagrangean density is given by 
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where 



B F (ip) = e(ip)~ 



(2) 



is the gauge kinetic function of a massless scalar field (p. 
This function acts as the effective dielectric permittiv- 
ity and can phenomenologically be taken as an arbitrary 
function of ip. 

Defining an auxiliary gauge field ajj = e A® and a new 
non-Abelian gauge field strength by 

1% = eF^ = 9 X ~ d "< + ^ abC «' (3) 
the covariant derivatives are written in the usual form 



D^ a = d^ a + e Q e abc a b l $ c , 

where e abc is the Levi-Civita tensor. 

We consider that the potential V($ a ) is given by 



(4) 



(5) 



where 7 > is the coupling of the scalar self-interaction, 
and v is the vacuum expectation value of the Higgs field. 
The Lagrangian density in cqn. Q is then invariant un- 
der SU(2) gauge transformations of the form 
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where A a is a generic isovector. This symmetry is broken 
down to U{1) because there is a non vanishing expecta- 
tion value of the Higgs field. Thus the two components 
of the vector field develop a mass Mw = veo, while the 
mass of the Higgs field Mh = v y/j. 

It is convenient to define the dimensionless ratio 



c 



M w e ' 



(8) 



and to rescale the radial coordinate by Mw, so that dis- 
tance is expressed in units of M w x . 

Varying the action with respect to the adjoint Higgs 
field <I> a ^ one gets 



Variation with respect to a° leads to 

a, [B F (<p)m =r 

with the current j° M defined as 

= e e a6c $ 6 (D"$ c ). 
Finally, variation with respect to <p gives 



in which f 2 = f"f atM - '. 
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We are now interested in static, spherically symmetric, 
magnetic monopole solutions. Therefore we make the 
"hedgehog" ansatz 



a?(r) = tiak 
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where x a are the Cartesian coordinates and r 2 
H(r) and W(r) are dimensionless radial functions which 
minimize the self-energy i.e., the mass of the monopole 
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where the coordinate r and the functions H and tp have 
been rescaled as 
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v ip. 



(17) 



It will prove useful to compute the energy density which 
is given by 
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with a prime meaning a derivative with respect to the di- 
mensionless coordinate r. Substituting the ansatz given 
in (fT5f> - <fT~5f> into equations [|§))- ljl"2l one gets 
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Defining 



<j> a 1 
f = _ f a j : 

|*| e |$| 3 



e a6c $ a (£> M $ 6 )(A / * c ) (22) 



and choosing the gauge where <!> a = 5 a3 |$|, i.e., with $ 
pointing in the same direction everywhere, one gets 



fnv = df, al - d v a 3 



(23) 



By writing a 3 = a M , one identifies (|23|l with the usual 
electromagnetic tensor which for constant ip satisfies the 
ordinary Maxwell equations everywhere except in the re- 
gion where H ~ 0. 
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III. NUMERICAL IMPLEMENTATION OF 
EQUATIONS OF MOTION 

In order to solve numerically the equations of motion, 
we have to reduce them to a set of first order equations. 
For that purpose we define the variables 



V = H 1 , U = W' , 77 = tp' 
Equations i|19l21fl can then be written as 
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FIG. 1: The numerical solution of the field H (r) as a function 
of distance, r, to the core of the t'Hooft-Polyakov monopole 
for several values of £. For £ = the Higgs field is massless 
and for £ — ♦ oo the Higgs field is frozen at its vacuum value 
except at the origin. 



which require at least six boundary conditions. 

At far distances from the core (r — > oo), the Higgs 
field H(r) falls off to its vacuum value, H = 1. Us- 
ing this boundary condition in eqn. I|19l) . one gets that 
W(r) must vanish far from the core. On the other hand, 
since the symmetry at the core is not broken, then the 
Higgs field vanishes and from regularity of the energy- 
momentum tensor, one can choose a gauge for which 
W = 1. The other two boundary conditions come from 
the normalization of the electric charge at the origin. At 
the core e = eo, which means that the gauge kinetic func- 
tion is equal to one. Finally, by using the Gauss law to 
solve the eqn. (|21() without sources of a variation other 
than the monopole, one gets that rj must vanish at r = 0. 
Therefore, we have at all six boundary conditions, four 
of them at origin and the other two far away from the 
core: 
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As the boundary conditions are at different points of 
the domain of the functions to be found, we use the re- 
laxation numerical method replacing the set of differen- 
tial equations by finite-difference equations on a grid of 
points that covers the whole range of the integration. 



A. t'Hooft-Polyakov Standard Solution 

In this section we consider as a first example the spe- 
cial case oi Bp = 1. This recovers the standard t'Hooft- 
Polyakov monopole problem [30l | in which there is no vari- 
ation of the fine-structure constant with ip(r) = constant. 



Let us first consider £ = 0, i.e., the case of a massless 
Higgs field. The energy of the monopole can be written 
as 
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where we have introduced a new scalar field P(r) as 

P{r) = H (1 - W 2 ). (32) 

Note that for H(r) and W(r) that solve the first-order 
differential equations 



H' = 

w = 

the energy is minimized to 

E = An [P(r -i 



1 - W 2 
-WH 



d) - P(r = 0)] 



(33) 
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This procedure was first pres ented by Bogomoln'yi (3jJ 
and generalized in j3£| [ifj, [41| . Taking into account the 
boundary conditions for H{r) and W(r), one gets that 
the minimum of energy is equal to 47r/eo. The solutions 



(36) 
(37) 



of equations (|33I34J) are 

H(r) = 

W(r) = 



sinh (r) 

1 1 
tanh (r) r 



which is the well known BPS (BogomoPnyi, Prasad, Som- 
merfield) solution [HIH]. In Fi § s - fflandJlwe plot H{r) 
and W(r) given by eqns. i|3l)|) and |J57J|. 

Let us now consider the opposite limit, £ — > oo. Fixing 
eo one gets in this limit that 7 — > 00, i.e., the Higgs 



FIG. 2: The numerical solution of the field W(r) as a function 
of distance, r, to the core of the t'Hooft-Polyakov monopole 
for three values of 




FIG. 3: The numerical solution of the energy of the t'Hooft- 
Polyakov E(£) as a function of the ratio of the Higgs to vector 
field mass, £. Note that the energy remains finite when £ — > 
oo. 



potential is much larger than the kinetic term forcing the 
Higgs field to be frozen at its vacuum value everywhere 
except at the origin. Then the only equation of motion 
is 



W" 



w 



W(l - W 2 ) 



(38) 



In Figs, n and [21 we plot H(r) and W(r) for £ — > oo. 

Of course, the former values of Q simplify very much 
the equations of motion. However, we have found nu- 
merically the set of solutions for a generic non vanishing 
finite C t . In fact, we achieved very good precision for Q 
ranging the interval 10 -4 < C < 10 3 - 

We also compute the mass of the monopole and ver- 
ified that it increases monotonically with increasing £, 
remaining finite when £ — > oo (see Fig. |3J). 



FIG. 4: The numerical solution of the scalar field ip = Ine 
as a function of distance, r, to the core of the monopole, in 
the Bekenstein model. Note that if ui < 0, e — > oo when 
r —* oo. On the other hand, if ui > then e — > when 
r — > oo. The dashed line represents the constant-a theory, 
which corresponds to the limit u> — > oo. 



Bekenstein Model 



Consider a gauge kinetic function, Bp(ip), given by: 



Bp{ip) = e 



(39) 



where w is a positive coupling constant. This recovers 
the original Bekenstein model|lg], 

Defining ip = tp/y/tjj and using eqn. l|3*9"|) the La- 
grangean density in can be written as 



C 



-2i)> 



(40) 



Note that despite the fact that the gauge kinetic function 
in eqn. I|39|l is only well defined for uo > 0, the model 
described by the Lagrangean density in eqn. I|40() allows 
for both positive and negative values of to. 

In this model equations 12(j|) and (|27|l are given by 



U' 



2aU 
2a 



W(W 2 
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e 

UJ 



-2ip 



e 2 ^'H 2 W (41) 
2X 2 



1 - W 



U 2 



(42) 



where we took a = tp'. First we note that in the limit 
uj — > oo one recovers the t'Hooft-Polyakov classical so- 
lution described in the previous subsection. In Fig. 0|we 
plot the numerical solution of the scalar field ip( r ) = lne 
for several values of u>. 

From Fig. 0|one concludes that if lu < then e diverges 
asymptotically away from the core of the monopole, and 
the energy density 



-2 i/i 
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C. Polynomial Gauge Kinetic Function 




FIG. 5: The numerical solution of the energy density p as a 
function r, the distance to the core of the monopole. The 
dashed- line is for uj — > oo. Note that for uj = 2.0, the en- 
ergy density is non negative, while for uj = —0.1 it becomes 
negative after some distance of the core. 




FIG. 6: The numerical solution of the Higgs field H(x) and 
the gauge field W(x) as a function of distance r, to the core of 
the monopole for. Note that the change in W(r) with respect 
to the t'Hooft-Polyakov solution (dashed-line) is more relevant 
than the change in H(r). 



H' 2 



(WH\ l 1 9n , 
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is no longer positive definite. However if u > then e 
vanishes asymptotically when r — > oo, and the energy 
density is in this case positive definite (See Fig. E}. We 
also note by observing Fig. 21 that in the large uj limit 
the curves for positive and negative uj are approximately 
symmetric approaching the dashed line which represents 
the constant-a model when to — > oo. 

Finally, in Fig. H3 we plot the numerical solution of the 
scalar field H (r) and the gauge field W{r) as a function 
of distance, r, to the core of the monopole. The dashed- 
line represents the constant-a solution and the solid line 
represents the Bekenstein one for uj = 2.0. We verified 
that even in the uj — > limit, the change in H(r) with 
respect to t'Hooft-Polyakov solution is still negligible. 



We now consider another class of gauge kinetic func- 
tions given by 



N 



B F (y) = l.0 + Y,PiV\ 

i=l 



(44) 



where Pi are dimensionless coupling constants and N is 
an integer. 

Note that by considering 



Pi 



(-2)* 



w 



i/2j\ 



(45) 



one recovers the Bekenstein coupling given in l|39ll . This 
relationship between the coupling constants Pi and uj 
has interesting consequences for the model given by eqn. 

First we verified that the behaviour of ip both for 
fli > and (3i < is similar to that of the Beken- 
stein model with uj > which is recovered in the limit 
of small | Pi | /large uj. 

Another feature that we verified is that if one takes 
Pi = one gets the t'Hooft-Polyakov limit, for any of 
Pi for i > 1. This means that there is a class of gauge 
kinetic functions for which the classical static solution is 
maintained despite the modifications to the model. 

Another property can be noted when one substitutes 
the gauge kinetic function in equation of motion for ip 
(|21) . One gets 

+ c 2 ff^AO/W*" 1 ) 



(46) 



with 



C 



,. 2 , 1 /(1-tF 2 ) 



(47) 



Since that C 2 > 0, when Pi — > — Pi one sees by eqn. 
(|46|l that (p(r) — > —cp(r). However, as Bp in eqn. I|44|l 
is kept invariant, H (r) or W(r) do not vary. 

Although we have found the set of solutions for several 
values of N, for simplicity we consider that N — 2 in eqn. 
that is 



B F = 1.0 + piuj + p 2 y 2 



(48) 



with two free parameters. We define the models 1 and 2 
as Pi = — 3, P2 = (linear coupling) and Pi — —2, p 2 = 
5 (quadratic coupling) respectively. 

In Fig. 0we plot the numerical solution of the scalar 
field ip(r) as a function of the radial coordinate. As we 
have shown earlier the model represents any model with 



FIG. 7: The numerical solution of the scalar field if) = Ine 
as a function of distance, r, to the core of monopole, for a 
polynomial gauge kinetic function. Models 0, 1 and 2 are 
defined by fix = (02 arbitrary), fix = —3, fi% — (linear 
coupling) and f3\ — — 2, (3% = 5 respectively. 



FIG. 10: The energy density as a function of the distance, r 
to the string core for models 0, 1 and 2. The dashed line rep- 
resents the constant-a model. Note that there is an increase 
of the energy density due to the contribution of the extra field 
ip in the case models 1 and 2. 
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FIG. 8: The numerical solution of the fields H(r) and W(r) as 
a function of distance, r, to the core of monopole, for models 
0, 1 and 2. Note that the change in W(r) with respect to the 
standard constant-a result is much more dramatic than the 
change in H(r). 



0.02 




FIG. 9: Plot of Li(r) = log(H,/H ) for the models 0, 1 and 
2. One clearly sees that even a small value of /3i leads to a 
different vortex solution from the standard t'Hooft-Polyakov 
one. 



[3i = and we have verified that the replacement (3\ —* 
—fix does not modify the solution for ip. 

One clearly sees in Fig. [S] that the change in W(r) 
with respect to the standard constant-a result is much 
more dramatic than the change in H(r), we define the 
function 

Li(r)=)og(j±Y (49) 

and plot in Fig. [HI the results for the different models. 

Note that even a small value of f3\ will lead to a mod- 
ification of the magnetic monopole solution with respect 
to the standard t'Hooft-Polyakov solution. 

We have also studied the behaviour of the energy den- 
sity in this model. In Fig. 1 101 one sees by comparing with 
the dashed line, which represents the constant-a model, 
that since the fine structure constant varies there is a 
new contribution due to the field ip, to the total energy 
of the topological defect. In fact the energy density of the 
monopole can be divided into two components: one that 
is localized inside the core of the monopole and other re- 
lated to the contribution of the kinetic term associated 
with the spatial variations of the fine structure constant. 

IV. CONSTRAINTS ON VARIATIONS OF a 

After having discussed in the previous section the nu- 
merical solutions for the varying-a monopole, our next 
investigating point is to find an overall limit of the spa- 
tial variations of the fine-structure constant on monopole 
networks. 

For simplicity we assume 

B F {tp) = l+/3 ip, (50) 

i.e., the gauge kinetic function is a linear function in ip 
that satisfies the spherically symmetric Poisson equation 
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given by 



V 2 <p=-0f. 



(51) 



Note that (3 is constrained by Equivalence Principle tests 
to be such that \(3\ < 1CT 3 G 1 ' 2 (see H2). 

Integrating ean. (|50() from the core up to r max , which 
represents a cosmological cut-off scale, one gets 



, dip 



W-f = /3I(r)M(r max ). 
dr 



(52) 



In ean. H52|l we used the mass of the monopole which is 
given by 



M(r) = 4tt / p(r')r u dr' 
Jo 



and 



J(r) 



M(r max ) J 



f 2 (r'y 2 dr' , 



(54) 



which is a slowly varying function of r outside the core 
always smaller than unity. Thus we can take I(r) ~ const 
and integrate eqn. to get 



f3IM(r max ) (I 1 



47T 



(55) 



where ro is a integration constant which could be identi- 
fied as the core radius. Using Bp(ip) = e~ 2 one gets 



1 



(3 2 IM(r n 



1 1 

r r 



(56) 



which means that the variation of the fine structure con- 
stant away from the monopole core is proportional to the 
gravitational potential induced by the monopolcs. 



where we have used a — q.q/ ' Bp(ip), with a — a{r — > oo) 
and «o = a ( r o)- 

B. Planck Monopoles 

Proceeding as above for the Planck scale symmetry 
with 



M{r 



10 19 Gev 



(59) 



one obtains an overall limit for the spatial variation of 
the fine-structure seeded by magnetic monopoles 



— < io- 7 , 

a 



(60) 



(53) 

which is still very small even for Planck scale monopoles. 



V. CONCLUSION 

In this paper we investigated static monopole solutions 
in the context of varying-a theories based on Bekenstein- 
type models. First we studied various models with 
constant a and reviewed the standard static t'Hooft- 
Polyakov magnetic monopole solution. Then we consid- 
ered models with varying-a and showed that despite of 
the existence of a class of models for which the t'Hooft- 
Polyakov standard solution is still valid, in general, our 
solutions depart from the former one. We showed that 
Equivalence Principle constraints impose tight limits on 
the variations of a induced by magnetic monopoles. This 
confirms the difficulty to generate significant large-scale 
spatial variation of the fine structure constant found in 
previous works, even in the most favorable case where 
these variations are seeded by magnetic monopoles. 



A. GUT Monopoles 

Let us estimate an overall limit for the spatial variation 
of a outside the core of GUT monopoles. In this context 
the mass of the monopole is of order of 

M( W ) ~ 10 16 GeV. (57) 

The variation of a from the core up to infinity is 
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